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We review the basic ideas of the dynamical mean field theory (DMFT) and some of the insights 
into the electronic structure of strongly correlated electrons obtained by this method in the context 
(3J[), of model Hamiltonians. We then discuss the perspectives for carrying out more realistic DMFT 

studies of strongly correlated electron systems and we compare it with existent methods, LDA and 
\ LDA+U. We stress the existence of new functionals for electronic structure calculations which allow 

us to treat situations where the single-particle description breaks down such as the vicinity of the 
(N ■ Mott transition. 



I. INTRODUCTION 



The last two decades have witnessed a revival in the study of strongly correlated electron systems. A large variety 
of transition metal compounds, rare earth and actinide based materials have been synthesized. Strong correlation 
effects are also seen in organic metals, and carbon based compounds such as Bucky balls and carbon nanotubes. 
These systems display a wide range of physical properties such as high-temperature superconductivity, heavy-fermion 
behavior, and colossal magnetoresistance to name a few ffl. 

Strong correlation effects are the result of competing interactions. They often produce at low temperatures several 
thermodynamic phases which are very close in free energy, resulting in complex phase diagrams. As a result of these 
competing tendencies, strongly-correlated electron systems are very sensitive to small changes in external parameters, 
O , i.e. pressure, temperature, composition, stress. This view is supported by a large body of experimental data as well 

as numerous controlled studies of various models of strongly-correlated electron systems ^, Q . 
,_H \ At the heart of the strong-correlation problem is the competition between localization and delocalization, i.e. 
J> ■ between the kinetic energy and the electron-electron interactions. When the overlap of the electrons among themselves 
is large, a wave-like description of the electron is natural and sufficient. Fermi-liquid theory explains why in a wide 
range of energies systems, such as alkali and noble metals, behave as weakly interacting fermions, i.e. they have a 
Fermi surface, linear specific heat and a constant magnetic susceptibility and charge compressibility. The one-electron 
spectra form quasi-particle and quasi-hole bands. The transport properties are well described by the Boltzmann 
^vq ■ theory applied to long lived quasi-particles, the approach that makes sense as long as kfl 1. Density functional 
theory (DFT) in the local density or generalized gradient approximations (LDA or GGA), is able to predict most 
"p-; ■ physical properties with remarkable accuracy. || 

When the electrons are very far apart, a real-space description becomes valid. A solid is viewed as a regular array 
of atoms where each atom binds an integer number of electrons. These atoms carry spin and orbital quantum numbers 
giving rise to a natural spin and orbital degeneracy. Transport occurs via activation, namely the creation of vacancies 
and doubly occupied sites. Atomic physics calculations together with perturbation theory around the atomic limit 
allows us to derive accurate spin-orbital Hamiltonians. The one-electron spectrum of the Mott insulators is composed 
of atomic excitations which are broadened to form bands that have no single-particle character, known as Hubbard 
• • ■ bands. In large number of compounds spin and orbital degrees of freedom order at low temperatures breaking spin 
rotation and spatial symmetries. However, when quantum fluctuations are strong enough to prevent the ordering, 
. possible new forms of quantum mechanical ground states may emerge Q. 

These two limits, well separated atoms, and well overlapping bands, are by now well understood and form the 
basis of the "standard model" of solid-state physics. One of the frontiers in strongly correlated electron physics 
problem is the description of the electronic structure of solids away from these limits. The challenge is to develop 
new concepts and new computational methods capable to describing situations where both itineracy and localization 
are simultaneously important. The "standard model" of solids breaks down in this situation, and strongly correlated 
electron systems have many anomalous properties, such as resistivities which far exceed the Ioffe-Regel-Mott limit 
PMott 1 ~ (e 2 /h)kj^, non-Drude-like optical conductivities, and spectral functions which are not well described by 
the band theory Q. 

To treat these systems one needs a technique which is able to describe Kohn-Sham bands and Hubbard bands on the 
same footing, and which is able to interpolate between the atomic and the band limit. Dynamical mean-field theory 
(DMFT)Q is the simplest approach satisfying these requirements. We introduce it in a very general formulation, the 
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cellular DMFT or C-DMFT ||, which is particularly well suited for electronic-structure calculations in section II . 

The goal of these lectures is to introduce recent DMFT developments to both the electronic-structure community 
as well as researchers interested in the many-body physics of correlated materials. For the electronic-structure 
community, DMFT is a promising technique for going beyond the LDA method. To illustrate the promise of the 
technique we describe in section III some of striking progress which has been achieved in the theory of the Mott 
transition by the use of DMFT at the level of model hamiltonians. 

The view of strongly-correlated electron systems that we describe in this introduction, stresses the need for incor- 
porating electronic-structure methods in treating strongly correlated electron systems. The low temperature physics 
of systems near localization-delocalization crossover, is non universal, system specific, and very sensitive to the lattice 
structure and orbital degeneracy which are unique to each compound. We believe that incorporating this information 
into the many-body treatment of these systems is a necessary first step before more general lessons about strong- 
correlation phenomena can be drawn. The extreme sensitivity of the materials properties to microscopic details has 
motivated us to realistic studies of correlated solids within DMFT. To put these efforts in an electronic-structure 
perspective, and to stress the qualitative difference between DMFT and other electronic structure methods, we review 
in sections IV and V the density functional method and the LDA + U method from an effective action point of view. 
In the following sections VI— XII we rely heavily on the effective-action formulation of dynamical mean-field theory 
for electronic-structure calculations || 1Q 



In section X we argue that many DMFT results obtained so far are in much better agreement with experiments 
than the corresponding results of LDA calculations We conclude in section XII with a brief introduction to the E- 
DMFT method, a complementary approach to C-DMFT which can take into account the longer range of the Coulomb 
interactions. This method, which we call the GWU method generalizes both the well established GW method 
and the DMFT but has not yet been implemented in a realistic framework. 



II. CELLULAR DMFT 

Reference^ reviews the pre 1995 work on the dynamical mean-field method and its various extensions. In this 
section we discuss a recently proposed || cellular version of the theory or C-DMFT. This formulation which is well 
suited for electronic-structure calculations since it is adapted to a non-orthogonal basis. This supercell (or cluster) 
DMFT remains close in spirit to the DMFT ideas, where the clusters have free (and not periodic) boundary conditions. 
Furthermore, the flexibility of this approach stresses the connection between the lattice many-body problem and self- 
consistent impurity models as in the single site dynamical mean-field theory The construction is carried out in 
complete analogy with the standard dynamical mean-field construction , but allows the use of a large class of basis 
sets. This frees us from the need to use sharp division of space into supercells. 

It has been proved || that the C-DMFT construction is manifestly causal, i.e. the self-energies that result from 
the solution of the cluster equations obey S(k, w) < 0, eliminating a priori one of the main difficulties encountered 
earlier in devising practical cluster schemes. 

It is useful to separate three essential elements of a general DMFT scheme: (a) Definition of the cluster degrees of 
freedom, which are represented by impurity degrees of freedom in a bath described by a Weiss field matrix function 
Go- The solution of the cluster embedded into a medium results in a cluster Green's function matrix and a cluster 
self-energy matrix, (b) The expression of the Weiss field in terms of the Green function or the self-energy of the 
cluster, i.e. the self-consistency condition of the cluster scheme, (c) The connection between the cluster self-energy 
and the self energy of the lattice problem. The impurity solver estimates the local correlations of the cluster, while 
the lattice self-energy is projected out using additional information, i.e. periodicity of the original lattice. 

Our construction applies to very general models for which lattice formulation naturally appears. It can be thought 
of as an extension of the band-structure formalism that takes into account the electron-electron interactions. The 
lattice hamiltonian, H[fi a , fj a ], (one example could be the well-known Hubbard hamiltonian) is expressed in terms 
of annihilation and creation operators /j CT and f\ a where i runs over the sites of a <i-dimensional infinite lattice 
i = (ix, . . . , id), the index a denotes an internal degree of freedom such as a spin index or a spin-orbital or band index 
if we consider an orbitally degenerate solid. 

(a) Selection of cluster variables: The first step in a mean-field approach to a physical problem, is a selection of a 
finite set of relevant variables. This is done by splitting the original lattice into clusters of size T[j=i Lj arranged on 
a superlattice with translation vectors R n . On this superlattice we choose wave functions \R n a) partially localized 
around R n with a = 1, . . . , N denoting an internal cluster index. The relation between the new wave functions, \R n a), 
and the old ones, \icr), is encoded in a transformation matrix, Sr^o^ct, such that \R n ot) = ^2 ia . N 17 )^^ 1 ^ a - Due to the 
translation symmetry of the lattices we have Sn na .ia = S aa (r(i) — R n ) where r(i) is the position of site i. The creation 
and annihilation operators of the new basis are related to the operators of the old basis by CR na = ^R n a,iafia- 
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and the operators that contain the "local" information that we want to focus our attention on are c a = C(R n =o)ai 
i.e. the operators of the cluster at the origin. We will refer to these operators as the cluster operators. Note that 
we do not require that the wave-function basis is orthogonal, and the nonorthogonality is summarized in an overlap 
matrix O™" = 0^(R m - R n ) = (R mf i\R n is). 

The next step is to express the hamiltonian in terms of the complete set of operators CR mPl . In terms of the new 
set of variables it has the form 



H — — 2_j t^v{Rm — Rn) c R mf i C R n v 

+ U^p.dR^C+^C+^CR^CR 



3 P 



(1) 



We stress again the generality of the method. Equation (Q) has the form one would obtain by writing the full 
hamiltonian of electrons in a solid in some tight -binding non-orthogonal basis. The hamiltonian is then split into 
three parts, H = H c + H c b + Hb where H c involves only the cluster operators and their adjoints, Hb contains CR nfJt 
with R n ^= only and plays the role of a "bath", and finally H c b contains both CR raM with R ^ and the cluster 
operators Cp (which have R n = 0). Physically H c b couples the cluster with its environment. A similar separation can 
be carried out at the level of the action, in the coherent state functional integral formulation of this problem, where 
the partition function and the correlation functions are represented as averages over Grassman variables, 



Z = I II Dc R na Dc Rna e 

R n a 



(2) 



where the action is given by 



f d A E 4. m pO™d T CR nV + H[4 m ^c Rn A 

Jo \R m pR n v ) 



s c + s, 



cb 



Sh 



(3) 



The effective action for the cluster degrees of freedom is obtained conceptually by integrating out all the variables 
cr„p with R n ^ in a functional integral to obtain an effective action for the cluster variables c M = CR m= op, i.e. 

Z 



U n 



(4) 



R m ^0,p 



Note that the knowledge of the exact S e ff allows us to calculate all local correlation functions involving cluster 
operators. As described in 0], this cavity construction if carried out exactly would generate terms of arbitrary high 
order in the cluster variables. Our approximation renormalizes the quadratic term, and neglects the renormalization 
of the quartic and the generation of higher order terms. Since the action S c b contains only boundary terms, the effects 
of these operators will decrease as the size of the cluster increases. Within these assumptions, the effective action is 
parameterized by Go !(II /(t — r'), the Weiss function of the cluster and has the form 



S eff - - 



f ! drdr'J24(r)G^l u (r^r')c^) 
J o pv 

I o\t x oIt 2 dr 3 dr 4 r Atiyps c+ (r x ) c+ (t 2 )Cs(t 4 )c p (t 3 ) 
Jo 



(5) 



where T^p^ = ^^-({O}). Using the effective action (||) one can calculate the Green functions of the cluster 



G c ,pu{T — t')[Gq 



-(T t Cp(t)c+(t'))[Go] and the cluster self energies 



= G 1 — G c , 



(G) 



(b) Self-consistency condition: The cluster algorithm is fully defined once a self consistency condition which indi- 
cates how Go should be obtained from E c and G c is defined. In the approach that we propose here the self-consistent 
equations become matrix equations expressing the Weiss field in terms of the cluster self-energy matrix E c . 



(7) 
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where O(k) is the Fourier transform of the overlap matrix, i(k), is the Fourier transform of the kinetic energy term of 
the hamiltonian in Eq. (|l|) and k is now a vector in the reduced Brillouin zone (reduced by the size of the cluster, Lj, in 
each direction). Equations (^) and (0) can be derived by scaling the hopping between the supercells as the square root 
of the coordination raised to the Manhattan distance between the supercells and generalizing the cavity construction 
of the DMFT ||] from scalar to matrix self energies. If the cluster is defined in real space and the self-energy matrices 
could be taken to be cyclic in the cluster indices so that the matrix equations could be diagonalized in a cluster 
momentum basis, Eq.(^) would reduce to the DCA equation ]50| . However, in the DMFT construction, the clusters 
have free and not periodic boundary conditions, and we treat a more complicated problem requiring additional matrix 
inversions. 

(c) Connection to the self-energy of the lattice: The self-consistent solution, G c and S c) of the cluster problem can 
be related to the correlation functions of the original lattice problem through the transformation matrix Sn m a,ia by 
the equation 

where S is the Fourier transform of the matrix S with respect to the original lattice indices i. Notice that "Siat.aa' is 
diagonal in momentum and will also be diagonal in the variable a if this variable is conserved. 

(d) Connection to impurity models: As within the single-site DMFT it is very convenient to view the cluster action 
as arising from a hamiltonian, 

Himp — ^ ] Ep^Cp C s + ^ ^ ^ pi/ps c p C v CpC^ 
pt; pvpt, 

+ e v a ij a v + ( v v,» a ij c » + h - c ) ■ ( 9 ) 

Here e^j is the dispersion of the auxiliary band and Vkj ;M are the hybridization matrix elements describing the effect 
of the medium on the impurity. When the band degrees of freedom are integrated out the effect of the medium is 
parameterized by a hybridization function, 



A^(iu; n )[e kj , Vkj] = ^2 , kj,/f _ kj '* . 



(10) 



The hybridization function is related to the Weiss-field function by expanding Eq.(|7|) in high frequencies: 

Ga 1 (iu n ) = iu n O-E-A(iu n ) (11) 

with O = [53k O k 1 ] 1 indicating that the impurity model has been written in a non-orthogonal local basis with an 
overlap matrix O 

Finally we note that one can view the C-DMFT as an alternative to the usual approach to treating finite systems 
by imposing on them periodic boundary conditions. Here we use a boundary condition which is physically equivalent 
to embedding the cluster in an infinite medium, which is determined self consistently using information computed 
within the cluster, and structural information on how the cluster is embedded in the infinite lattice. 

III. QUALITATIVE INSIGHTS FROM DMFT APPLIED TO MODEL HAMILTONIANS 

1. DMFT phase diagrams, Frustration, Complexity and Universality 

The low-temperature phase diagram of simple hamiltonians treated within DMFT has several distinct phases, and 
is fairly complex. Even the simplest, bare-bones hamiltonian (one-band Hubbard model with partial frustration) has 
at least a metallic antiferromagnetic phase and a paramagnetic insulating phase in addition to a paramagnetic metal 
phase and the antiferromagnetic insulating phase. The understanding of this model required several years of research 
efforts by several groups [£3| |24| ]2j| |6l| |Q ||| pi] . In this section we highlight some of the insights obtained, to 
motivate the need for extending the DMFT method to incorporate realistic aspects of the electronic structure. 

The phase diagram) 14 shown in Fig [l] shares some similarities with the phase diagram of Vanadium oxide. This 



observation leads Rozenberg et. al.[h4 to suggest several optical experiments which have confirmed some qualitative 
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predictions of DMFT. These and other successful predictions of DMFT applied to simple model of correlated electrons 
are described in this section to motivate the applications of this technique in a more realistic setting. 

It is important to emphasize, however, that the main lesson drawn from the qualitative similarity in the low- 
temperature region between the DMFT phase diagram of one of the simplest models of correlated electron systems 
and that of some real oxides, is the ability of the DMFT method to capture multiplicity of possible ordered states. 
The detailed nature of these phases and the character of the transitions between them depend on many details of 
the hamiltonian describing the specific crystal structure and chemistry of the compound. To approach this problem 
realistic versions of DMFT have been constructed and are being developed jl5|, |l6| . 

Strong dependence of the low-temperature phases and of the low-temperature physical properties of each material 
on its crystal structure and chemical composition should be contrasted with the remarkable degree of universality 
that is predicted to occur at higher temperatures. All that is required to produce the high-temperature features of 
the DMFT phase diagram is a large degree of magnetic frustration to suppress the long-range order and to allow for 
a localized phase with a large entropy content. In systems without magnetic frustration, the onset of magnetism or 
other forms of order preempts us from accessing this strongly correlated regime. The origin of the magnetic frustration 
is crucial for understanding the low-temperature part of the phase diagram, with its myriad of ordered phases, but 
is rather irrelevant in the high-temperature regime, where thermal fluctuations average all various configurations 
leading to a more universal description which is captured by a relatively local approach such as DMFT in its single 
site or in its clusters versions. In systems such as titanates and vanadium oxides, the origin of frustration arises from 
the orbital degeneracy which is unique to those materials. In nickel selenide sulfide mixtures, the crystal structure is 
such that a sizeable ring exchange term competes with the nearest-neighbor superexchange interaction resulting in a 
reduced Neel temperature. Still, these systems display very similar phenomena around the Mott transition endpoint. 

Contrast between highly-universal behavior at high temperature and the dependence of low-temperature properties 
on additional parameters in the Hamiltonian, was discussed [p"l[ in connection with the comparison of the physical 
properties of the vanadium oxide and the nickel selenide sulfide mixtures. The phase diagram of the two-dimensional 
organic compound k BEDTTF ]l7| , where the frustration originates in its underlying chiral triangular lattice of 
dimers, strengthen the validity of this point of view. Indeed many of the high-temperature physical properties of this 
have been accounted for by the DMFT studies of McKenzie and Merino [lq] . 

To summarize, since magnetic frustration and competition of kinetic and interaction energy are all that is required 
for obtaining the high-temperature part of the "canonical" phase diagram of a correlated electron system. At low 
temperatures between two values of U , U c \ and U C 2 two dynamical mean field solutions are possible. The transition 
between the localized and extended regime as a function of j- takes place via a first-order transition |6^| Jl4j , this is 
faithfully reproduced by the simplest model containing these ingredients treated within DMFT. 
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FIG. 1: Schematic phase diagram of partially frustrated Hubbard model from ref the inset illustrates the behaviour of the 
resistivity above but near the Mott endpoint 



The phase diagram|14| displays two crossover lines. The dotted line in fig 1 is a coherence incoherence crossover (i.e. 
the continuation of the l7 c2 line where metallicity is lost). The shaded area is a continuation of the U c i line, where the 
temperature becomes comparable with the gap. Both were observed in the V2O3 and NiSeS system jl9[ |2Cj. Further 
justification for this point of view, and a refined description of the localization derealization transition around the 
Mott transition endpoint was achieved by the development of a Landau like description plj, |22fl . 



6 



2. Coherent and Incoherent Spectra 



Mapping onto the Anderson impurity model offers an intuitive picture of both metallic states and Mott insulating 
states. A correlated metal is described locally as the Anderson impurity model in a metallic bath: The Kondo effect 
gives rise to strongly renormalized quasiparticles when the interactions are strong, and to a broad band when the 
interactions are weak. The Mott insulator is locally described as the Anderson impurity model in an insulating bath. 
The charge degrees of freedom are gapped, but the spin degrees of freedom are not quenched; they dominate the 
low-energy physics. When there is one electron per site, the Mott transition takes place as one goes from the first 
regime to the second by increasing the strength of the interaction U ||| . 

A sketch of the evolution of the spectral function —ImG(iVL + iS) of the half-filled Hubbard model is described in 
figure (§) . For interactions U close but smaller than the critical U C 2 the one-electron spectral function of the Hubbard 
model in the strongly correlated metallic region contains both atomic features (i.e. Hubbard bands) and quasiparticle 
features in its spectra |]l2| . This may be understood intuitively from the Anderson- Yuval-Hamman path-integral 
representation. In the regime of strong correlations, paths that are nearly constant in imaginary time, as well as 
those that fluctuate strongly, have substantial weight in the path integral. The former give rise to the Hubbard 
bands, while the latter ones are responsible for the low-frequency Kondo resonance. Two features in these spectra 
are surprpizing. First, the narrow central peak before the Mott transition, resulting from quasiparticle states formed 
in the background of the coherent Kondo tunneling of the local spin fluctuations. Second, atomic physics leaves a 
signature on the one-particle spectral properties in the form of well-formed Hubbard bands at higher frequency, in 
the strongly-correlated metallic state. 

As the transition at zero temperature is approached, there is a substantial transfer of spectral weight from the low- 
lying quasiparticles to the Hubbard bands. The Mott transition at zero temperature takes place at a critical value of 
U, denoted by U C 2 where the integrated spectral weight at low frequency vanishes, as shown in figure pi This results 
in a Mott transition point where the quasiparticle mass diverges, but a discontinuous gap opens in the quasiparticle 
spectra^] . These results are in agreement with the early work of Fujimori et al. [ p6[ , who arrived at essentially the 
same picture on the basis of experiment. It is worth remarking that the spectral function in the strongly correlated 
metallic region is better regarded as composed of three components: Hubbard bands centered at U/2, a low-energy 
quasiparticle peak with a height of order unity and a total intensity proportional to (U C 2 — U) distributed over an 
energy range U C 2 — U , and an incoherent background connecting the high energy to the low-energy region. The last 
feature ensures that there is no real gap between the quasiparticle features and the Hubbard bands, as long as one 
stays in the metallic regime. 
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FIG. 3: pdc{T) around the coherence incoherence crossover near the finite temperature Mott endpoint. U/D = 2.1,2.3,2.5 
(bottom to top), obtained with the IPT method from ref |l4|. 



Anomalous resistivities 



Figure || describes the anomalous resistivities near these crossover regions. Notice the anomalously large metallic 
resistivity which is typical of many oxides 0. While the curves in this figure far exceed the Ioffe-Regel limit (using 
estimates of kf from T=0 calculations) there is no violation of any physical principle. At low temperatures, a k space 
based Fermi-liquid-theory description works but in this regime the resistivity is low (below the Ioffe-Regel limit). 
Above certain temperature the resistivity exceeds the Ioffe-Regel limit but then quasiparticle description becomes 
inadequate. There is no breakdown or singularities in our formalism, the spectral functions remain smooth, (above 
the Mott-transition endpoint), only the physical picture changes. At high temperatures we have an incoherent regime 
to which the Ioffe-Regel criteria does not apply, because there are no long-lived excitations with well-defined crystal 
momentum in the spectra. The electron is strongly scattered off orbitals and spin fluctuations, and is better described 
in real space. In this regime, there is no simple description in terms of k space elementary excitations, but one can 
construct a simple description and perform quantitative calculations if one adopts the spectral function as a basic 
object in terms of which one formulates the theory. 

Only the anomalously large magnitude of the resistivity (which follows from the Green function which has branch 
cuts rather than well-defined poles), is universal as can be seen by comparison of the detailed temperature dependence 
at half filling (as in fig ||) and away from half filling as in figures U and |[ The temperature dependence of the transport 
in the high temperature incoherent regime depends on whether the system is at integer filling or doped, as can be 
shown numerically Jl|| and analytically 27 in the example of the doped Mott insulator. The low-temperature and the 
high-temperature anomalously large resistivities also occur in strongly coupled electron-phonon systems, as discussed 
earlier lEil. 



4- Anomalous Transfer of Spectral Weight 



Another manifestation of the same physics is the anomalous transfer of spectral weight which is observed in the 
one-electron and in the optical spectra of correlated systems as parameters such as doping or pressure are varied. This 
surprising aspect of strong-correlation physics was noted and emphasized by many authors |£9j. Transfer of spectral 
weight can also take place as a function of temperature. For example the "kinetic energy" which appears in the 
low-energy optical sum rule can have sizeable temperature dep endence, an effect that was discovered experimentally 
by [[30) and explained theoretically by DMFT calculations |5l) 

Once more, thinking about this problem in terms of well-defined quasiparticles is not useful. It is more fruitful 
to formulate the problem in terms of spectral functions describing on the same footing coherent and incoherent 
excitations. The relative weights of these components in the spectra evolves smoothly with temperature and leads 
to sizeable variations in the integrated optical intensity. The evolution of the spectral function near the temperature 
driven Mott transition is shown in figure ^|. 

The physics of strongly correlated materials in a wide range of parameters cannot be described in terms of quasipar- 
ticle excitations. Recent advances in the theory of the Mott transition highlight the fundamental role of the spectral 
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FIG. 4: pdc(T) in units of zha/e 2 , vs T (in units of D) for different values of orbital degeneracy N for a fixed doping 5 — .1 
obtained with the NCA method which is valid at high temperatures, from ref 
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FIG. 5: PdJT) vs T in units of zha/e 2 , vs T (in units of D) obtained with the IPT method, for different dopings at ^ = 2.8 



from ref 27 



function. The functionals discussed in sections VII, VIII, and XII carry over this ideology to realistic electronic- 
structure calculations. 

IV. DENSITY FUNCTIONAL THEORY AND LDA 

Density functional theory in the Kohn-Sham formulation is one of the basic tools for studying weakly-interacting 
electronic systems as is widely used by the electronic-structure community. We review it briefly here using an 
effective action point of view in order to highlight the similarities and differences with the DMFT methods which 
will be presented in the same language. The approach in this context was introduced by Fukuda |3^, J&j, and 
it amounts to a simple Legendre transformation. One introduces the generating functional of the connected Green 
functions in the presence of an arbitrary potential. 

Consider the partition function Z (or equivalently the free energy W) of a system of electrons moving in a crystal 
potential V ex t{x) and interacting via Coulomb interactions V, in the presence of an external source J coupled to the 
electron density 



Z = exp[-W[J}} 



D[ip^]exp[-S- / dxJ{x)^ a \x)i> a (x)] 



(12) 
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(T-T c )/T c =-1 .8% 

(T-T c )/T c =-0.9% 

(T-T c )/T c = 0% 

(T-T c )/T c =+0.9% 

(T-T c )/T c =+1 .8% 




FIG. 6: Evolution of the spectral function as function of temperature (bottom to top), near the finite temperature from Ref. 
p^ . The inset is the spectral function at the second order Mott endpoint. 



S= I ' dxY,^i{x)[d T -^ + V ext (x)]ip r7 (x) 

a 

+ f dxdx'i>+(x)i,+ (x')U c (x - x')i> a ,{x')^ a (x) 



(13) 



In (|l^), x = (x, r) denote space-imaginary time coordinates, V ex t is the crystal potential created by the ions and 
U c (x — x') = |x — x'\~ 1 5(t — t') is the Coulomb interaction. 
The density functional T(p) is defined by 



T[p] = W[J] - / J(x)p(x) 



(14) 

where one eliminates the source in favor of the density. The minimum of the density functional gives the true 
density and the total energy of the solid. 

To construct approximations to the functional T it is very useful to introduce the Kohn-Sham potential, Vks, 
which is defined as the potential such that added to the non-interacting kinetic energy produces the given density in 
a reference system of non-interacting particles . i.e. 

The exact functional can be viewed as a functional of two variables 

T(p, V KS ) = -Tj2*r log[iu n + V 2 /2 - V KS ] - J V KS (r)p(r)dr+ (16) 

\ J ^-r'l dTdr ' + J V e*t(r)p(v)dr+E xc [p] 

T(p) is obtained by substituting Vks(p) obtained by solving Eq. ( |l5| ) (which makes (|l^ ) stationary) onto T(p, Vks)- 
E xc [p] is the exchange-correlation energy which is a functional of the density and not of the external potential. 
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Extremizing (]l^) with respect to p gives 



1 A ■ <tr]J ' 1 = / + Vxc{r)[p] + K:rt(r)W (17) 



where 14 c (r) is the exchange-correlation potential obtained as 

5E XI 



c . , - ^xc(r) (18) 

Since -E^clp] * s n °t known explicitly approximations are needed. The LDA assumes 

E xc [p] = J e xc [p{r)]p(r)dr (19) 

with e xc [p(r)] being the energy density of the uniform electron gas, a function which is easily parametrized. 

For simplicity we restrict ourselves to zero-temperature and interpret the Fermi functions as the Heaviside step 
functions. Eq. (|l5|) can be rewritten as the eigenvalue problem 

[-V 2 /2 + Vfcs(r)hMr) = e kj -^(r) (20) 
p(r)=^/(e kj )Vk J -(r)Vk j (r) (21) 

kj 

and Vks is given as an explicit function of the density. In practice one frequently uses the analytical formulae fl35|. 
The idea here is to fit a functional form to a Quantum-Monte-Carlo (QMC) calculation, and we will return to this 
idea when we discuss the IPT. 

Then, the total energy of the crystal is given as 



E to t = Y,f^j) e ^ + \ J P \r-^\ dTdT ' + J V ^t(r)p{r)dr+ J e xc \p{v)]p{v)dv + E dc 

kj 



(22) 



E dc = - V KS (v)p{r)dr (23) 



where 



simply subtract the interaction energy from the Kohn-Sham eigenvalues which is explicitly included in the Hartree 
and exchange-correlation term to avoid double counting. 

Density functional is an exact approach as long as the invertibility condition expressing the potential in terms of the 
density is satisfied. However in strongly correlated situations, the total energy is not very sensitive to the potential 
since the electrons are localized due to the interactions themselves, and the lack of sensitivity of the functional to 
the density, hampered the design of good approximations to the exact functional in this regime. Furthermore, at the 
Mott transition the invertibility condition may not be satisfied. 

The effective action approach exhibits that the functional T is a Legendre transform of the exact generating func- 
tional W[J] and introduces Kohn-Sham field in a perturbative expansion of the Legendre transform in e 2 at fixed 
density. 

The Kohn-Sham parametrization of the density in terms of Vks is extremely useful because it expresses uniquely 
p(r) in terms of the Kohn-Sham orbitals ipkj(r). Truncations of DFT are necessary for practical implementations. 
There are two different philosophies in truncating Eqs.(^)— (22), one is via the introduction of pseudopotentials which 



we will not describe here. All-electron methods simply introduce a finite basis set Xa( r ) ancl expand 

lMr)=£xiS( r )4? (24) 

a 

keeping a finite set of a. Notice that this truncation immediately restricts the active part of the multiplicative operator 
associated with the Kohn-Sham potential to have a form 

V = Y,\X«) V <*P(X*\ (25) 
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Of course, one can add to this contributions from the set which is orthogonal to the minimal basis set |Xq) without 
changing the truncated density. The requirement of locality presumably determines the Kohn-Sham potential and 
the component of ( |25| ) outside the space ofy^(r) uniquely. 

Linear muffin-tin orbitals (LMTO's) |36| are an optimal minimal basis set. For a known Kohn-Sham potential 
this construction can be done once and for all. However, since Vks depends on the density, the basis is 
adapted iteratively to the self-consistent solution. The same observation applies to our new implementation of the 
LDA+DMFT approach. 

In principle, density functional theory is an exact theory as long as the invertibility condition expressing the potential 
in terms of the density is satisfied. In practice, the usefulness of this approach is due to the existence of very successful 
approximations to the exchange-correlation potential, by the LDA or the GGA. In principle the Kohn-Sham orbitals 
is a technical device for generating the total energy, however in practice they are used as a first step in perturbative 
calculations of the one-electron Green function, as in the GW method. The LDA method is very successful in many 
materials for which the standard model of solids works. However, in correlated electron system this is not always 
the case. Our view, is that this situation cannot be remedied by using more complicated exchange and correlation 
functionals in density functional theory. As discussed in the previous section the spectra of strongly correlated electron 
systems have both Hubbard bands (which have no analog in one-electron theory) and correlated quasiparticle bands 
in the one-electron spectra. DMFT is the simplest approach which is based on this simple physical picture. 

The extension to spin-density functional, in its non-collinear version is straightforward. The functional T should 
now be considered as a functional of both the charge density p{r) and magnetization density m(r) The latter is a 
vector quantity. It can be slightly non-collinear due to spin-orbit coupling effects. The functional of the total energy, 
Eq.(p2|) is rewritten by taking into account the fact that the Kohn-Sham field now consists of both the potential 
VKs( r ) and magnetic field Bks(^) 



Etot = ^/(eki)e kj + ^ P ^^ drdr' + J V ext (r)p{r)dr + J B ext (r)m(r)dr (26) 
e xc [/j(r),|m(r)|] x p(r)dr+ J f xc [p(r), |m(r)|] x |m(r)|dr + E dc 



where we have used the local spin density approximation (LSDA) expression for exchange-correlation energy by 
assuming that the properties of electron gas depend on the absolute value of the magnetization: |m(r)|. 
The double-counting term is 

E dc = - J V KS (r)p(r)dr- J B^ s (r)m(r)dr (27) 

with 

fiE 

B KS (r) = B ext (r) + ( 28 ) 
om(r) 

m(r) = ^fkj{ifkj\s\'ipkj) S pin (29) 

(we average over spin degrees of freedom only), and spinor wave functions f/'kj satisfy to the Pauli-like Kohn-Sham 
matrix equation: 



[(-V 2 + V KS (r))I + B KS (r)s + #8]^ = e kj ^ kj 

where / is the unity 2*2 matrix; s is a spin operator which is just the Pauli matrix divided by two; £ determines the 
strength of spin-orbit coupling and in practice is determined [fj7| by radial derivative of the I = component of the 
Kohn-Sham potential inside an atomic sphere: 

ew - m 

cr ar 



V. LDA+U METHOD 



We now turn to a description of the LDA+U method^]. We will deal with the spin unrestricted formulation and 
for illustration purposes assume that relativistic effects are small. This allows us to choose a quantization axis along 
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some direction, say z, since the total energy is now invariant with respect to its orientation. Instead of considering 
p(r) and full vector m(r) we can deal with p(r) and only m z (r) or, equivalently with spin-up and -down densities 
p°(r) = [p(r) + am z (r)]/2 , a = ±1. 

The approach requires an introduction of a set of localized orbitals </> a (r — R) which are used to build an "occupancy 
spin density matrix" 

Kb = E fM I ^,»<Mr)dr f ^- CT (r')^(r')dr' (31) 

This intuitively represents the "correlated part of the electron density" as long as we associate our projectors <p a with 
correlated electrons . The total energy now is represented as a functional of the spin densities P a M and of n° b . In 
complete analogy with Eq. (|l^) , one introduces a Lagrange multipliers matrix A^ b such to enforce (pi] ) and to expand 
the LDA+U functional as 

T LDA+U [nl b ,\l b ,V£ s ,p°] = -T^Vlog^n + V 2 /2 _ ycr s - ^ Kbta^Wb (*')) ~ (32) 



(33) 



a J a ab J \ \ 

E Model [n a ] -E^ odel [n a ] 
where we have added a contribution from the Coulomb energy in the shell of correlated electrons 

E Mo de i K] = I £ £ U abcd n° ab n-J + \ E E " -W<f><* (34) 



<T abed cr afccd 



Since, partially this energy is already taken into account in LDA, we have to subtract a double-counting part denoted 
by Ef c odel [n a ]. 

Minimization of this functional gives rise to 



iuJ n \ 

with 



[iw„ + V 2 /2 - VZ S - E KbMrWbir'T 1 



ab 



r \ e -»0+ = ^ /(ek . £r) |^. ff(r) |2 (35) 



[-V 2 /2 + ^ S + E W.W^M]^. = ekjabve (36) 
where V£ s is given by equation similar to Eq.([T7|) and 

Kb = E U abcdn^ + }XUgbcd - Jabcd)n a cd ds — (37) 

cd cd ab 

Several remarks are in order. 

(i) The LDA+U functional and the LDA+U equations are defined once a set of projectors and a matrix of interac- 
tions U a bcd is prescribed. Formally, the matrices U and J have the following definitions: 

U a bcd = (ac\-\bd) = I (t>\{T)r c {r')v G {T-T')<t>b{v)(t> d {r')drdT' (38) 

' J a 

Jabcd = (ac\-\db) = f <p a ^r)<f>* c (r')vc(r-r')<fi d (r)Mr')drdr' (39) 

where the Coulomb interaction i>c(r — r') has to take into account the effects of screening by conduction electrons. 
In practice, one can express these matrices via a set of Slater integrals which, for example, for d-electrons are 
parametrized by three constants F^°\F^ 2 \ and F^. 

An important question is to discuss the double counting term E M od el [n a ] which has been occurred when we added 



additional Coulomb interaction to the functional. It was proposed |39| that the form for E^ I c odel [n a \ is 
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Ef c odel = \pn(n - 1) - 1 J[flT(flT - 1) + - !)]• 



where 



U 

J = U 



2/(2 ,Vi)Bh> & )-h>-)) 

ab 



(40) 

(41) 
(42) 



and where n* 7 = ^ a n^ a , and n = rzJ + n - ' - . This generates the correction to the potential in the form: 

Kb = E UabajnZ + Y,V<*<* ~ J >*cd)<d ~U{n-\) + J{ff - \) ( 43 1 

cd cd 



The question now arises as to whether the double-counting correction should subtract self-interaction effects or 
not. In principle, if the density functional contains this spurious term, the same should be taken into account in the 
double -counting expression. Judging by the experience that the LDA total energy is essentially free of self-interaction 
(the total energy of the hydrogen atom is, for example, very close to -1 Ry, while the Kohn-Sham eigenvalue is only 
about -0.5 Ry), the construction E d VI c odel is made so as to be free of the self-interaction. However this statement 
cannot be considered seriously in general, and alternative forms of the double counting term includiing the effects of 
self-interaction have been used pql 

(hi) If one uses the Eqs.(pTj), (|37j) and ( |3l| ) to eliminate V£s, n^ b , and \° b as functions of p a , and backsubstitutes 
that into (|32"|), one obtains a functional of the density alone. 

(iv) The literature is ambiguous on whether Tlda+u could perhaps be viewed as an approximation to a functional 
of a projected density matrix , or, as indicated in (iii), simply as a different density functional where in addition 
to dividing the density into spin up and spin down as in LSDA, one introduce a correlated component (i.e. what is 
parametrized by n a b) of the density and an uncorrelated one. 

(v) A different point of view is to introduce a " correlated part of the one particle density matrix" 



a (r)(^ + (r)V(r'))4(r')rfrdr' 



and to consider a functional of this quantity and of the total density by effective action methods and view Eq. 
as an approximation to this exact functional. However, the interpretation of 



(44) 



V 2 /2 



V KS 



a b 



(or 1 



(45) 



as a density matrix is not consistent with existence of interactions This is because Eq. (jMj) describes a density 
matrix which has eigenvalues which are ones and zeros, and this characterizes a non-interacting density matrix. The 
density matrix of an interacting system has eigenvalues which are less than one. Therefore, Eq.(]3l|) can not represent 
an interacting density matrix. Because of these difficulties we have suggested to interpret the LDA+U method as a 
static limit of the more powerful DMFT method which we describe later. The static limit of DMFT is going to be 
most accurate as more symmetries (spin, orbital) are broken. Removing local degeneracies by spontaneous symmetry 
breaking is the simplest way of minimizing the energy (reducing the correlations) [{fH . 
(vi) It was argued Efl] that the Green function 



V 2 /2- 



ab 



\ ab MrWb(r')Y 



(46) 



can be viewed as a limiting case of the GW approximation |5J but this is again not clear since an interacting Green's 
function has poles with residues less than one, and this is not the case in expression (|46| ) except for the uncorrelated 
situation when Hartree Fock theory becomes exact. 

(vii) Relativistic effects which are important for such applications such as magnetic anisotropy calculations, can 
be considered. We have described the extended DFT in the previous section, and here we only discuss the LDA+U 
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corrections. If spin-orbit coupling is taken into account, the occupancy matrix becomes non-diagonal with respect to 
spin index: 



kj J 



^(r)0 o (r)dr / <(r')^(r')dr' (47) 



The correction to the functional has the form similar to Eq(34) and it is given by 



zp Model r„ a a' 1 ^ \ A j j ^aa „ — a— a , ^ \ A /r r t \ atr ^aa ^ \ ^ j a—a^—aa ( a o\ 

E t n J = 2 ^ U "bcdn ab n cd I - 2, (£4bcd - Jabcd)n ab n cd - - ^ Jabcdn ab n cd (48) 

abcda abcda abcda 

which can be worked out by considering a Hartree-Fock average of the original expression for the Coulomb interaction 

1 e 2 

- 2j2J<ao-&o-'| — \cada')c+ 7 c& T ,c da >c ar (49) 

(Ttr' abed 

The correction to the potential takes the same form as Eq.(fl3|) when a = <t', i.e., 

Kb = E U -bcdn-J- a + ^iUabcd - JabcdXl ~U{n-\) + J(n™ - ~) (50) 

cd cd 

and for the off-diagonal elements it is given by 

Kr = -J2 J ^ n cd aa (51) 

cd 

To make it more physically transparent we can introduce magnetic moments at the given shell by 



s 



/ j aa 



M Kb (52) 



where fi runs over x,y,z for Cartesian coordinates, or over,— 1,0, +1 (z, ±) for spherical coordinates. Relativistic 
correction to the LDA+U energy can be written in physically transparent form 

^ JabcdK^n-™ =\Y1 m aV Jabcdm^ + \Y1 ^b J abcdm ( ^ ] (53) 

abcda abed abed 

In principle one can assume further generalization of exchange matrix J a bcd to be anisotropic, i.e depending on fi/j,': 

J abcd- 

To summarize, since the density uniquely defines the Kohn-Sham orbitals, and they in turn, determine the oc- 
cupancy matrix of the correlated orbitals [once a choice of correlated orbital in Eq. (147])] is made, we still have a 
functional of the density alone. However it is useful to proceed by analogy with Eq. (|32|), and think of the LDA + 
U functional as a functional of n a V£ s and A CT , whose minimum gives better approximations to the ground-state 
energy in strongly correlated situations. Allowing the functional to depend on the projection of the Kohn-Sham 
energies onto a given orbital, allows the possibility of orbitally ordered states. This is a major advance over LDA in 
situations where this orbital order is present. As recognized many years ago, this is a very efficient way of gaining 
energy in correlated situations, and is realized in a wide variety of systems. 

Those are the formal difficulties of the LDA+U method. From a practical point of view, despite the great successes 
of the LDA+U theory in predicting materials properties of correlated solids (for a review, see Ref. ^) there are 
obvious problems of this approach when applied to metals or to systems where the orbital symmetries are not broken. 
The most noticeable is that it only describes spectra composed of Hubbard bands. We have argued in the previous 
sections that a correct treatment of the electronic structure of strongly correlated electron systems has to treat both 
Hubbard bands and quasiparticle bands on the same footing. Another problem occurs in the paramagnetic phase 
of Mott insulators, in the absence of any broken symmetry the LDA + U method reduces to the LDA, and the gap 
collapses. In systems like NiO where the gap is of the order of eV, but the Neel temperature is a few hundred Kelvin, 
it is unphysical to assume that the gap and the magnetic ordering are related. For this reason the LDA+U predicts 
magnetic order in cases that it is not observed, as, e.g., in the case of Pu 0|. 
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VI. MODEL HAMILTONIANS AND FIRST PRINCIPLES CALCULATIONS 



It is useful to think of a program for performing realistic electronic structure calculations for correlated materials 
in the light of the qualitative discussion of the earlier section. The hamiltonian describing electrons at short distances 
is known and easily written down. This is the formal starting point of all-electron first-principles calculations. So, 
the theory of everything is summarized in the action S of a system of electrons with Coulomb interactions between 
them, described in by the action of eq. (|l3|). 

This action ignores relativistic effects, which are important for the structural properties of heavy elements. Rela- 
tivistic corrections which introduce spin-orbit interaction have the form Hps — Hso + H$oo + H$s 




HsOO = — r ij X Pi\ s i + 2s j) 



a 2 m-^ 1 (nj-Sj), 

i<j y y 

In addition this action should be supplemented by the electron-phonon interaction which can be not negligible in 
correlated materials. They have been investigated using DMFT within model calculations j5?J , but realistic studies 
are in its infancy and will not be considered here. 

Introduction of a complete tight binding basis <^>z(x) allows us to rewrite the action in the form 

S = J dT[cJO a ^d T c p + H[cjc p ]] (54) 

with the Hamiltonian containing an infinite number of bands, 

ff = - E E E 4 m + h - c -i + E E E^'V c ^ c ^ c ^4.' ( 55 ) 

ij Im a iji'j' klmn aa' 

the new operators c are related to the continuum operators ifi by c ia = J V , o-(x)0i(x — Rj). k,l,m,n denote band 
indices, aa' the spin indices and ij denote the lattice sites. O describes the overlap matrix if this basis is not 
orthogonal. 

However, to describe the physics at a lower energy scale one would like to eliminate the degrees of freedom which 
have energies much larger than that scale, and derive an effective Hamiltonian which is more transparent and contains 
only the relevant or active degrees of freedom. The effective Hamiltonian at that scale, is the model Hamiltonian 
which is usually written down by the solid-state physicist on physical grounds. 

The explicit construction can be written down as a Wilsonian elimination of the irrelevant high energy degrees of 
freedom. Formally, one divides the set of operators in the path integral in ch describing the "high energy " orbitals, 
and cl the low-energy orbitals that one would like to consider explicitly. 

exp(-S e ff[c L c L ]) = J Dc H Dc H exp- S[c H ,c L ,c L ,c H ] (56) 

In the electronic-structure program this procedure is called downfolding and is carried out at the level of Kohn-Sham 
orbitals. The transformation (|56| ) generates retarded interactions of arbitrarily high order. However there are good 
reasons, why an approximation in which only quartic terms are kept, and its frequency dependence ignored is an 



excellent approximation. This results in models of the form (55) except that the interactions and the hoppings are 
screened relative to their bare values. 

If one keeps the constant part of the effective action, which contains the free energy of the high-energy orbitals 
which have been eliminated this procedure in principle contains all the information which is needed to carry out total 
energy calculations. Nothing is lost except for the higher-order interactions and the retardation which is unlikely to 
be very important. The hamiltonians obtained in this way are the model hamiltonians used in many-body theory. 
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Finally, the elimination in Eq ( J56| ) cannot be performed even within the approximations cited above. However any 
technique which can be used to treat Hamiltonians approximately, can also be used to performed the elimination (|5(f). 
In particular carrying the local density approximation has been applied to obtain the U's. It would be interesting to 
reconsider the validity of these approximations from the point of view of an effective action. 

As the number of orbitals and the relevant energy scale is reduced a renormalization group (RG) flow in the space 
of all hamiltonians is defined. Different initial conditions at short distances describe different substance, or materials 
at different pressures, lattice spacings, dopant concentration etc. 

If one starts with conditions that correspond to weakly correlated systems ( e.g. atomic numbers involving s or p 
electrons, high densities etc. ) the RG flows are relatively simple and converge at low energies to reach simple fixed 
points describing band metals or insulators. 

On the other hand when we start from more correlated situations (e.g. open shells, containing relatively localized d 
or / electrons, lower densities), the RG trajectories are diverging from one another, reflecting the diversity of phases 
nearby. This situation calls for quantitative methods for realistic modelling of the material in question. One of the 
most serious difficulties in carrying out the Wilson RG program described above, is the continuous change in the form 
of the effective Hamiltonian from scale to scale. A typical example is the formation of a heavy fermion liquid state at 
a coherence energy scale. At high energies the effective Hamiltonian contains atomic configurations and conduction 
electrons, at low frequencies only heavy quasiparticles are the relevant degrees of freedom. In spite of these difficulties, 
an R.G analysis taking account some quantum chemistry in the initial conditions has been carried out in the local 
approximation aided by developments in DMFT ]4^| 

While following the R.G. flows down to very low temperatures and predicting physical properties in the most 
strongly correlated situations may prove to be very difficult, there are many reasons to believe that C-DMFT with 
small sizes will be accurate in a wide range of interesting situations (not too close to phase transitions, not too low 
temperatures). 

The process of eliminating degrees of freedom with the approximations described above gives us a physically rigorous 
way of thinking about effective Hamiltonians with effective parameters which are screened by the degrees of freedom 
to be eliminated. Since we neglect retardation, and terms of higher order than six the effective hamiltonian would 
have the form 



% = -EEE 4 m + h.c] + E EE v^ n ,4 a k 4^ aC } a , + e (57) 

ij Ira o iji'j' klmn aa' 

It should be regarded as the effective hamiltonian that one can use to treat the relevant degrees of freedom. If the 
dependence of Eq on the nuclear coordinates are kept, it can be used to obtained the total energy. If the interaction 
matrix turns out to be short ranged or has a simple form, this effective Hamiltonian could be identified with the 
Hubbard or the Anderson hamiltonians which have been treated in the literature. 

To conclude this section it is worth clarifying several terms that are used in electronic-structure literature, the 
light of the previous discussion, for pedagogical reasons. The first point is the meaning of ab initio or first-principles 
calculations. These imply that no empirically adjus table parameters are needed in order to predict physical properties 
of compounds, except the structure and the charges of atoms. First principles do not mean exact or accurate or 
computational inexpensive. If the effective hamiltonian is derived, i.e. if the functional integral is performed by a 
set of well-defined approximations, and the consequent hamiltonian (^7|) is solved for its total energy keeping track 
of the constant Eq we have a first-principles method. In practice, the derivation of the effective hamiltonian or its 
solution may be inaccurate or impractical in which case the ab initio method is not very useful. We make this remark 
because H e ff has the form of a "model Hamiltonian" and very often a dichotomy between model hamiltonians 
and first-principles calculations is made. What makes a model calculation semiempirical is the lack of coherent 
derivation of the form of the "model hamiltonian" and the parameters entering there. A second point has to do with 
the elimination of degrees of freedom and the meaning of all-electron calculation, indicating that all the electronic 
degrees of freedom are taken into account. All electronic structure methods make some elimination of degrees of 
freedom, it is most evident in pseudopotcntial methods, but it is also true in other methods. To make a problem 
in the continuum computationally feasible, a discretization is necessary. The elimination of unnecessary degrees of 
freedom greatly facilitates the accuracy of the consequent discretization. Introduction of a finite basis set in electronic- 
structure calculations is a form of discretization, the important question is its accuracy. The LMTO's have proved to 
be extraordinarily accurate to discretize the Kohn-Sham hamiltonian in the relevant region of energy. In DMFT not 
only the one-body term of the Hamiltonian is truncated, but also the interaction terms are, this truncation is most 
accurate in ultralocalized basis. The C-DMFT and the E-DMFT are many-body approaches which lend themselves 
to practical truncations in non-orthogonal LMTO's. 
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VII. DYNAMICAL MEAN-FIELD THEORY 

^From a conception point of view, the construction of dynamical mean-field functional constitutes a radical de- 
parture from Kohn-Sham based DFT. The dynamical mean-field equations as we will see are the equations for a 
continuous distribution of spectral weight. Attention has been shifted away from the well-defined Kohn-Sham quasi- 
particles (poles in the Green function) to continuous distribution of spectra (which appear as branch cuts) and give 
rise for example to the Hubbard bands. 

Another central difference connected to this point is the fact that the local spectral function can now be identified 
with the observable one-electron spectrum. This is very different from the Kohn-Sham quasiparticles which cannot 
be identified rigorously with the one electron spectra. While DMFT is computationally more demanding than DFT, 
it is formulated in terms of observables and gives more information than DFT. It can be formulated in an effective 
action point of view in complete analogy with DFT. 

Our starting point is a multiband Hubbard Hamiltonian, i j denote lattice sites while 

a /3 are spin -orbital indices: 

H = E E l t4^ + E E < 15 ^ 4 ( 58 ) 

ij a0 ij a/3fS 

which is obtained from the procedure discussed in the previous section. The effective action construction of DMFT 
parallels that given in other sections. A source Ji a f3{T,T') is introduced in the partition function 

e -0wm = J dc + dce -s~E ia , i j ia ^y)4 a (T) Cifi (T')dTdT> (59) 

W generates the local Green function 

SW i \ 

j- j^ = -(c ia (T)c+(T>)) = A aP (ry) (60) 

The effective action is T(A) = W[j(A)] - J (A) A. The minimum 

gives rise to correct spectra and total energy. Notice that the definition of the local Green function depends on the 
basis of orbitals chosen, this is similar to the LDA+U method. However, unlike the LDA+U approach, we are dealing 
with the functionals of Green functions which have a meaning independently of the Kohn-Sham representation of the 
density. 

Notice that the definition of the local Green function requires the choice of a basis. However, unlike the LDA+U 
method, the choice of basis is used to perform a Legendre transfomation with respect to a well-defined object, the 
local Green function, rather than with respect to a part of the density, as in the LDA+U method which lacks a clear 
physical significance. 

The starting point (|5|) can be regarded as a model hamiltonian, but, as argued in the previous section, if the 
constant parts are kept and (|57]) is carefully derived, this is effectively equivalent to the full hamiltonian in the 
relevant energy range. 

The dynamical mean-field approximation to the functional T, i.e. Tdmft, can be written in two alternative forms 
depending on whether we stress that it is a truncation of the exact functional when expanding T in powers of the 
hopping (atomic expansion) or in powers of the interaction (expansion around the band limit). To write this functional 
it is useful to define the quantity x(iw n , A). The physical meaning of this expression is parallel to the meaning of the 
Kohn-Sham potential: it is the function that one needs to add to the free hamiltonian in order to obtain a desired 
spectral function, as described in table I. 



A a P{iw n ) = Yfon ~ *(k) - eo - X(«"n, A)]~l (62) 
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TABLE I: Parallel between the different approaches, indicating the physical quantity which has to be extremized, and the field 
which is introduced to impose a constraint (Kohn Sham field). 

Method Physical Quantity Constraining Field 

LDA plf) V KS (r) 

DMFT (band limit) A a/3 (iuj) Xaf)(iu) 

DMFT (atomic limit) A a/3 (iuj) A a @(iu)) 

Baym-Kadanoff G a /s(k, iu)) Y, a p(k,iu>) 

LDA+U p(r),n ab V K s(r),X a b 

LDA+DMFT (band limit) p{r), A a0 (iu) V K s(r),xMiu) 

LDA+DMFT (atomic limit) p(r), A a0 (iu) V K s(r), A afj (iu) 



This quantity is important to set up an expansion of the functional in powers of the interaction following the inversion 
method. This emphasizes DMFT as a partial summation of the expansion around the band limit. Notice that if the 
exact self-energy of the problem is k independent, then \ coincides with the self-energy, this statement is parallel 
to the observation within DFT: if the self-energy of a model is k and frequency independent then the self-energy 
coincides with the Kohn-Sham potential. The DMFT functional is given by 



T DMFT {A, x) = -J2J2 tr l °s(^n - t(k) - e - X) - tr X A (^n) + $[A] 



(63) 



where $[A] is a sum of all local graphs (on a single site i) constructed with V^ lS as a vertex and A as a propagator 
which are two-particle irreducible. The diagrammatic rules for the exact functional are more complicated and were 
discussed in ref . [ ^) . 

Eq. ( |62| ) determining \ — x(A) appears as a saddle point of the functional (63) and should be backsubstituted to 
obtain T dm ft {A), the DMFT approximation to T(A). 

The second derivation emphasizes the expansion around the atomic limit, where the starting point is a dressed 
atom. It starts by writing the hamiltonian (p57|) into two parts: H = H + Hi with H = ^2h [i] , where 



i in 



khnn 



klmn +k 4-1 m n 
liiiaa' ia 1 ia ia ia' 



(64) 



and Hi = H — H Q and carries out the inversion method in powers of A-ffi (A is a coupling constant to be set to unity 
at the end). Here the indices k, 1, m,n are orbital indices. 

The zeroth-order term requires the introduction of the Kohn-Sham field, Jo = J^fj, which couples to the local 
lattice Green's function A [m The effective action involving H and Jo nas the form of an impurity model where the 
source lays the role of the hybridization therefore we will switch notation and sometimes use A(r, r') for this source. 
This action is highly non-trivial due to the presence of the impurity-interaction term in ( |64|) . The lowest-order term 
in the inversion method then requires us to solve an impurity model which expresses the source A(r, r') in terms of 
the impurity model Green functions of a generalized Anderson impurity model defined by 



S at [A] = J drdr' 



aim 



d 



S (T- T ')^7 + A /m (r-r') 



C(r') + / drH 



(65) 



The impurity Green function is given by A?£ (r, r') 



SW at 



7T = - (cia(T)c+ a (T')) = A lm (T,T') 



SJlm(T,T') '/,„ 

where the expectation value of any operator O is given by 

/ Dc+ l Dc™0 exp — 5 [ A] ] 



(O)j 



and 



/ Dc+ l Dc™cxp-S[A 
W at [A] = -log J dc+d 



(66) 



(67) 



S at [c + c] 
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describes an atom or a set of atoms in the unit cell embedded into the medium., a labels spin, and orbital including 
position in the unit cell. A is the Weiss field of the mean-field theory. It is simply the Kohn-Sham field (with respect 
to the expansion around atomic limit) which is defined by the equation expressing that A gives rise to the local Green 
function A (see Eq. (|6^) where the general Kohn-Sham source was identified with the hybridization of the Anderson 
impurity model). 

The functional corresponding to (|6^) is now given by 

T [A] = W at [A[A}} - A[A]A (68) 
= Tr log A — A~lA + $ [A] 

with the A~t = iuj n — e where e = ^ k e(k) (e(k) is the Fourier transform of the hoppings iy, W at the free energy 
A = and $ is the sum of all one particle irreducible diagrams constructed with the local vertex Van and A. 

Now the functional T[A] can in principle be constructed if all terms in the expansion in A could be summed up. 
The DMFT approximation is obtained by summing all diagrams with a topology of a cactus, describing multiple 
excursion away from a given site, leading to the expression 

T DMFT {A, X ,A) = W ot [A]-£V(AA)-£VlogA- ( 69 ) 

JF^2^2 tr lo g(* w n - *( k ) - £ o - X) - ^2 tr(x - iu n + e )A 

It is useful to check the equation obtained by differentiating with respect to A which results in the CPA condition: 

iuj n - e - A = A^ 1 + x{iu n ) 

The tr runs over the spin and orbital indices, N s is the number of unit cells 

One can rewrite this functional by expressing both A and x m terms of the Weiss field Jo (which in this formalism 
plays the role of the Kohn Sham field of the effective action formalism when an expansion around the atomic limit is 
carried out, as opposed to the expansion around the band limit considered by Aliev and Fernando). 



r[A[A], X [A]] = Wq[A] - Tr log A - £ Tr log[tw„ - e(k) - X ] - ^[x - iw„ + e + A] A (70) 

k iuj n i^n 

A second important point is that both the truncation of the expansion around the atomic limit or the truncation of 
the graphs in powers of V, indicate that Tdmft is going to be a poor approximation to T(A) when the interactions 
V are highly non-local. This lead to the introduction of a hybrid method where the light electrons are treated by 
LDA which contain screening effects while the heavy electrons are treated by the dynamical mean-field methods. We 
will review this approach in the next section. 

Notice however that an exte nsion of the DMFT formalism allow us to bypass the introduction of LDA altogether 



as we will show in section XII 



VIII. LDA+DMFT 

We now turn to the LDA+DMFT method |l5|, which has recently been implemented in a self consistent way 
p6| . On one hand, this approach can be viewed as a natural evolution of the LDA+U method to eliminate some 
difficulties discussed at the end of Section V. On the other hand, it can be viewed as a way to improve the DMFT 
approach, so as to bring in more microscopic details to an approach that had been extremely successful at the level 
of model Hamiltonians as described in section III. 

We derive the equations following the effective action point of view |^| .To facilitate the comparaison between the 
approaches discussed in the earlier sections we have tabulated the central quantities which have to be minimized, and 
the fields which are introduced to impose a constraint in the effective action method |B2| . As in the LDA+U method 
one introduces a set of correlated orbitals </> Q (r — R). One defines an exact functional of the total density p(r) and of 
the local spectral function of the correlated orbitals discussed before: 

A ab (iw,R) = - J 0l{r-R) (^(r,iu)i}+{r',iLu)) <f> b (r' - R)drdr' = - (c a (ioj,R)c+(ico,R)) (71) 

where indexes a, b refer exclusively to the heavy orbitals. 
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We now introduce sources for the density and the local spectral function of the heavy orbitals, l(r), and j ab (R., t, t') 

e -PW[l,j] = J dc + dce- S -I 'WpWdr-E^fi/ Oah{R,T,T')c+ R (T)c bR (T')dTdT' (J^ 

SW *) (73) 



<H(r 

SW 



Sjab{R, r, t') 



(c Ra (r)c+ b (r')) = A ab ( T ,r') (74) 



Then, the functional of both density and the spectral function is constructed by Legendre transform. This is an exact 
functional of the density and the local Green function, T(p,A), which gives the energy at the stationary point, and 
for which in principle, a pcrturbative construction can be carried either around the atomic limit or around the band 
limit following the inversion method. 

However, based on the evidence that even at the level of simple model hamiltonians DMFT describes accurately 
the properties of many systems and contains physics which is not captured by any other approach, and based on the 
remarkable success of LDA in treating weakly correlated electron systems, a useful approximation suggests itself, in 
a form of a LDA+DMFT approximation. 

The functional implementation corresponding to this approximation is given by F lda+dmft (ft, Vks^i A) which 
has the form 

T LDA+DMFT {p,V KS X,A) = -T^trlog[ioj n + Vy2-V KS -X ab (iuj n )<t> a (r)<l>* b {r')}- (75) 

J V K s(r)p(r)dr - ^ab(iu n )A ba (iu n ) + 

iuj n ab 

J V ext {v)p{v)dv + \ J E^l dvdv' + E^ A [p] + 
J2®[A ab {j,iu)] + $ DC 

3 

The field \ of the previous section is now denoted by S. Together with Vks they are the fields needed to constrain 
the density and the local spectral function of their correlated orbitals to their given values p and A. 4> is the sum of 
all two-particle irreducible graphs constructed with the local part of the interaction and $£>c is taken to have the 
same form as in LDA+U method, i.e. a simple local Hartree-Fock form with n ab — A ab (iu})e lUn<) As we argued 
before in a fixed tight-binding basis, —V 2 + Vks reduces to H TB (\t) and the functional ^ lda+dm ft for a fixed 
density and truncated to a finite basis set takes a form identical to the DMFT functional Eq. ( |70| ) in the section VII. 
Its minimization leads to the set of equations with Kohn-Sham potential defined by Eq. (|l7]) and 

S« 6 (^n)= *f , +ef fe C (76) 
dA ab (iuj n ) 

which identifies E as a self-energy of a generalized Anderson impurity model in a bath characterized by a matrix of 
levels 

^ = ^ C + E^ B (k) (77) 

k 

and a hybridization function A ab (iuj n ) obeying a self-consistency condition 



iuiO ab - e° ab - A ab (iuj n ) = E afc (iw„) + ^(iuj n O - e° - f(k) - E(iw n )) 1 

. k 

Finally, minimizing Eq. (|7ji) with respect to Vks indicates that p(r) should be computed as 



-l 

a, 



(78) 
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p(r) = T]T(r 



o n + V 2 /2 -Vks-Y, ZabiiUn^airMir'T 1 



r ) e 



which as indicated before, when truncating in a fixed set of orbitals becomes 



(79) 



(80) 



The solution of Eqs.(|76[)-(p0[) is carried out in double iterational loop described in Fig. 7. It is important to note 
that the choice of basis determines the quality of the truncation in a given energy range. 




FormV xc [P] 
Solve impurity model to get E(i'co) 



I 



Synthesize LMTO matrices H(k),0(k) 
Energy shifts of correlated orbitals, 
Add self-energy to form Green function 

G 1 (ico)=(ico+|Li)0(k)-H(k)- £(/oo)-E 



I 



Computer density p(r), local Green 
function and Weiss field: 
PW= ^n^apXa(r)G a ^iCO)Xp(r) 

G (ico)=[i:[(ico-Rt)0(k)-H(k)- L(ico)-E T] 
k + E(ico) 



i j 



FIG 7. Schematic description of the implementation of the LDA+DMFT loop carried out in r ef[TJ 

The basis can be gradually refined so as to obtain more accurate solutions in certain energy range. In principle this 
improvement is done by changing the linearization energies, and the experience from density functional implementa- 
tions could be carried over to the DMFT case. Finally notice that the rational for the counter term described in Eq. 
( f77| ) originates in comparison of spectra with the results of LDA+U calculations and deserves further investigations. 
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A different possibility which is computationally very practical, is to adopt the philosophy of the generalized IPT, 
introduced earlier|^| and to fix the double counting term by requesting the fulfillment of the Luttinger theorem, i.e. 
to chose 



-E a6 (0) 



(81) 



A version of this counter term, has been recently applied to iron with encouraging results |65 The ambiguity in 
the choice of double counting correction may turn out to be hard to resolve since one is substracting terms from a non 
diagrammatically controlled approximation. The ideas described in section XII, which results in a full Green function 
treatment of the many- body problem is perhaps the most promising route towards a non empirical resolution of this 
ambiguity. 



IX. TECHNIQUES FOR SOLVING REALISTIC DMFT EQUATIONS 



In practice the solution of the DMFT equations is more involved than the solution of the KS equations, which now 
appear as static analogs. Even for the simplest model hamiltonians, some care and judicious choice is required to 
obtain correct results. There are two central aspects of DMFT. One is the self-consistency condition, Eq. (f78|). This 
step is trivial for model calculations but becomes time consuming when realistic bands are introduced. It is done 
using tetrahedron methods. Second, one has to solve the impurity model. Fortunately, we can now rely on several 
years of experience to device reasonable approximations. First, at sufficiently high temperature QMC calculations are 
possible. Notice however the poor scaling with orbital degeneracy. The approach that we have advocated to carry out 
this step is very similar in spirit to the analytic parametrizations of V xc in LDA. One uses different approximations 
to the self-energy of the impurity model viewed as a functional of A(iw) in different regions of frequency. Eq]. 



X. SUCCESSES 



The early studies on V2O3 and NiSeS which were qualitative in nature, indicated the great potential of the DMFT 
approach. Another early contribution was a series of papers which display semiquantitative agreement on the 30% 
level on a series of bulk physical properties for the LaSrTiO system. The specific heat was considered |f52f , the 
magnetic susceptibility was computed, the optical conductivity was calculated EJ|, the Hall effect was obtained [g3j. 
Finally the thermoelectric power was focused[Q. On those fillings, the difference between the one band and the 
multiband situation for low-energy properties was estimated not to be very large for low energy properties given 
the uncertainties in the determination of the interaction parameters. Given the fact that only very simple tight- 
binding parametrizations were used in those works, and the fact that a large number of experiments were fit with the 
same value of parameters one should regard the agreement as satisfactory. The photoemission spectroscopy of this 
compound, as well as others, are not completely consistent with the bulk data, and several authors have argued that 
disorder and modelling of the specific surface environment is required to improve the agreement with experiment J55| . 

Given excellent agreement at the level of model hamiltonian calculations it was natural to incorporate the full band 
structure into DMFT to perform calculations of spectra which would incorporate more high energy physics, such as 
additional bands. This was first implemented jlM ff5[ 

In all these systems the LDA+DMFT approach gives results which are in better agreement with experiments than 
simple LDA calculations. We have recently argued argued that correlation effects are important for subtle effects such 
as the magnetic anisotropy. The easy axis of iron and Ni were obtained correctly [p6|. Lichtenstein and collaborators 
have obtained an excellent theoretical description of the spectra of iron and nickel ||58[ . 

A major recent development is the full self-consistent implementation of the DMFT+LDA approach. This devel- 
opment sheds a light on the mysterious properties of S Plutonium, and has resulted in new physical picture of the 
a phase of this material^). A second recent development is the development of the algoritms to perform transport 
coefficients which allowed the evaluation of thermoelectrical properties of correlated materials J60[. 



XL CONNECTION WITH LANDAU FUNCTIONAL 



It is useful to consider explicitly the limiting case, of a single band model on a lattice, and to evaluate the functional 
( |69| ) in the simple case of the semicircular density of states with half bandwidth D, a toy model which has been 
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studied intensively over the past few years. In this case Eq. ( |62] ) which expresses x,or z = iuo — e — xin terms of A, 
becomes A = -^\z — \/ z 2 — Z? 2 ]which is easily inverted to give z — t 2 A + J 4" 1 with 2t = D. 

The energy functional can be evaluated explicitly since l°g[ z — *fe] = [ z — V z 2 — D 2 } + Log[ ^ z+ ^ z ^ D ^ ] + iand 
reduces to 



\Y imp [A]+Tj2(zA-2A(iu>)A- 1) = W imp [A]+T^(tA - -^{iuj)) 2 - ^ Y 



where 



-f3Wim P =Logl df + dfe-^ <U-L, a ,aJ „i«^w*rW ( 82 ) 

and L Ioc [/+, /] = ///^[^ + + Uf+rff+fi. 

The stationary condition which relates A to A is such that the second term vanishes at stationarity, so as far as the 
stationary points are concerned this is equivalent to the Landau functional pjj which was considered in the context 
of the theory of the Mott transition. 



F LG [A] = -T £ + F mp [A] (83) 

UJ 

This functional can be understood by analogy to the Hubbard-Stratonovich construction of the Ising model free 
energy. 

0F La [h]=p£j-log[ch[20h]] (84) 

The first term in Eqs. (|8^) and (H) represents the energetic cost of forming the Weiss field, while the second terms 
are the energy gain of the local entity (spin or electron in the classical quantum case respectively) due to the presence 
of the Weiss field. 

It is important that the Landau functional is not just constructed by Legendre transformations, but it contains 
additional source terms which vanish at the saddle point. This flexibility in extending the Baym-Kadanoff construction 
to obtain functionals with more desirable properties such as (|3|) by adding sources were discussed in Ref . flltl . 



XII. EXTENDED DMFT, GW AND GWU 

In this final section, we discuss how one can avoid the introduction of LDA in our realistic DMFT framework, while 
still retaining an accurate treatment of the Coulomb interactions which should allow an accurate evaluation of the 
total energy. This is based on extended dynamical-mean-field ideas, which have been developed in a series of works 
p9[ |68[ . This development is particularly important to avoid the ambiguities connected with the substraction of the 
double -counting correction, discussed in section VIII. 

The starting point is the Hamiltonian, 

h = E E C c+ *»^ + E E rQ/35 Vc^c^ + E E E ^tK^^c^ m 

ij a{3 i aftyS ij ab a0^8 

The form of the Coulomb interaction is not completely general, and in particular it omits pair hopping and correlated 
hopping terms of the kind discussed by Hirsch. The corresponding action is given by 

5 = E/ 4«J^S% + / Hdr (86) 

which can be immediately transformed into an action which resembles that of an electron-phonon problem by in- 
troducing Hubbard-Stratonovic fields for decoupling the non-local part of the Coulomb interactions. The physical 
meaning of the fields i<fi a (j) are the electric potential created by the multipole of the charges in all the cells different 
fromj. Notice that B a ^(i, i) = by construction since Vij = 0. The Bose propagator in equation |S7| simply reproduces 
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to the action S = S e + S p + S ep 



the effects of the Coulomb interaction. After these transformations we can apply the extended dynamical approach 

S p = IJ2[ I B ab(ri, r'j)r(Tt)4> b (r'j)dTdr' (87) 



2 

ij 



Se = EE r^+cfrv* + ct^o^f - £ ctt^ (88) 

i af3~fS ij 

S^'EE rm a aP ctc lf3 dr (89) 

i aa/3 J ° 

The Baym-Kadanoff functional for this problem was derived in ref : 

r[0,G,n] = S[4>] + igJ2 Kp<t> a (t)Gp*(i, i) + TrlogG - Tr[G^ - G^]G (90) 

i 

--Triogn + -Trfs- 1 - n _1 ]n + $[G,U] (91) 

= ~ / dxdy^B-^x - y)Mv) (92) 

G, IT are the full fermion and boson Green function, and $[G, II] is the sum of all 2-particle irreducible diagrams 
constructed using the vertex (|89| ) with particle lines G, II. 

The extended DMFT equations are obtained by retaining only the local terms i.e. $(G,II) sa ^edmft 
=Yln &[Gii, Tin]. The importance of this truncation, is that its solution can be formulated as a dynamical mean- 
field theory, that is one can introduce an impurity model as in ref, the only difference is that now a Weiss field for 
the electrons and for the bosons has to be introduced. 



G(i,i,icj n ) = [G^-Hr^M) (93) 



Tl{i,i,iu n ) = [B- 1 + 2^]- 1 (i,i) 

a a/3 

Our case of long range Coulomb interactions corresponds to TTo{q, iio n ) — V(q). The one-band case was considered 
in ref ^[ It produces qualitatively new effects turning the Mott transition from second order to first order. A full 
numerical solution of this problem was only recently carried out for a model Hamiltonian | p7| , in a study of a Fermi 
Bose system. 

E-DMFT treats both electrons and the collective excitations (spin and charge fluctuations ) on a similar footing. 
The fermion (boson) propagators G and II are expressed in terms of self energies Hp(iu; n ) £.e(ia; n ) which are taken 
to be momentum independent, G~ 1 (iuj n ,q) — iu> n — t q ~Y 1 p(iuj n ),Tl{q,iuj n ) = — £b + B(q,iui n ).Th.e self-energies as 
well as other local quantities are computed from the local action: 

Sloe = /drdr'^ C +(r)G -;(r-rOc ff (rO-EI^( T ) n ^( T - T >^ T ')+ (95) 

<7 a 

fdr^ r a05l 4 a 4 c llCl& + z £ [ 1/T rm a aa p4 a ^dr (96) 

J q/3 7< 5 a/3 



and the parameters of the local action are determined by solving the E-DMFT self consistency conditions 
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Iio\^n) = [J2 v jtt rxr \ -L n( J ' + £ B {H , G }^„) 



s 



-l 

- * g - s F {n , Go} 1 



Go V») = E — i ^ fTT ?7T l + E ^{ n "' G }(*^) (97) 



The advantage of this formulation is that it is naturally combined with additional diagrams so that it resembles a 
natural extension of the GW method. Hence we propose the approximation $gwu — ^emdft + is the lowest 

order nonlocal correction (which is likely to be much smaller that the ordinary Hubbard term). 

$1 M)lEEE GtfGl s AfA2 s D% (98) 

i^j ab afiyS 

If indeed the local approximation is a good starting point, there should be not much difference between a full self 



consistent treatement or a perturbative lowest order treatment of the self energy derived from Eq. (98) 



XIII. OUTLOOK 

It is important to stress in this work devoted to methods for treating strongly correlated electron systems, the fruitful 
interplay between many-body theory of model hamiltonians and the field of electronic-structure calculations. The 
hard work in developing the dynamical mean-field technique has paid off handsomely in insights into real materials, 
and has already resulted in a major advance in electronic-structure theory, i.e. in our ability to predict physical 
properties of materials starting from first principles. Furthermore the technology developed for the solutions of model 
hamiltonians has been effectively transferred into the field of electronic structure. Simple calculations which gave 
the first spectra of the Hubbard model, using mapping onto impurity models |T^| , have been implemented with full 
realism using the tetrahedron metho d |]i"5| . Simple functionals which were developed to analyze the simplest possible 
localization-delocalization transition [fLO^ [21| , have been extended and implemented to yield total energies of solids JXg| . 
Interdisciplinary exchanges are a two-way street, and many of the cluster DMFT ideas which originate in advances 
in one-electron theory, hold great promise for the use in model hamiltonians. There is no doubt that the study of 
simple systems has been essential in our quest to control and understand complex materials 

The DMFT techniques will continue to be used in a much wider range of realistic problems, where correlation 
effects are prominent. In this context |6^| none of the basic ideas of the local approximation, requires a translation 
invariant. In fact DMFT has been applied to strongly disordered system which exhibit Anderson localization and to 
surfaces |34|. There are many physical reasons to believe that correlation effects in alloys and interfaces are much 
stronger than in the bulk of periodic solids DMFT can also be applied to finite systems such as complex molecules or 
molecular clusters. 

The next step in the development of DMFT as a realistic electronic-structure method is its implementation in 
a molecular-dynamic calculation, to predict structures without any a priori information. A second important and 
largely unexplored direction is the study of systems far from equilibrium. There have been substantial advances 
in understanding single impurity Anderson model far from equilibrium in the context of quantum dots, and this 
understanding will certainly allow major advances in the treatment of solids when combined with the standard 
DMFT self-consistency condition. 
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